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Abstract — In this paper we study the fast adiabatic-like ex- 
pansion of a one-dimensional Bose-Einstein condensate (BEC) 
confined in a harmonic potential, using the theory of time-optimal 
control. We find that under reasonable assumptions suggested 
by the experimental setup, the minimum-time expansion occurs 
when the frequency of the potential changes in a bang-bang 
form between the permitted values. We calculate the necessary 
expansion time and show that it scales logarithmically with large 
values of the expansion factor. This work is expected to find 
applications in areas where the efficient manipulations of BEC 
is of utmost importance. As an example we present the field of 
atom interferometry with BEC, where the wavelike properties of 
atoms are used to perform interference experiments that measure 
with unprecedented precision small shifts induced by phenomena 
like rotation, acceleration, and gravity gradients. 

Index Terms — Quantum control, time-optimal control, Bose- 
Einstein condensate 



I. INTRODUCTION 

A Bose-Einstein condensate (BEC) is the state of matter 
emerging in a dilute gas of weakly interacting bosonic atoms 
confined in an external potential when cooled to temperatures 
very near to absolute zero. More specifically, below a critical 
temperature a large fraction of the bosons occupies the lowest 
quantum state of the external potential, and the quantum 
effects become apparent on a macroscopic scale. This peculiar 
state of matter was first predicted by Einstein, when he 
extended the statistics of light quanta (photons), developed 
by Bose JT|, to material particles 0. 

Since its first experimental demonstration in 1995 0, 10], 
BEC has become a workhorse for atomic physics experiments. 
One premier example is the use of BEC for precision mea- 
surements in the context of atom interferometry [5 1, [6 1, where 
the wavelike properties of atoms are exploited to perform 
interference experiments that measure small shifts induced by 
phenomena like rotation, acceleration, and gravity gradients 
Q. For this kind of applications, the necessity to control and 
manipulate the BEC is ubiquitous. 

Among the various necessary control steps, a very important 
task is the ability to expand the BEC without exciting higher 
states, an undesirable effect that generates friction and heating 
|8|. The conventional way to achieve this goal is to change adi- 
abatically the potential that confines the BEC. The drawback 
of this approach is the long necessary times which may render 
it impractical. Recently, a method to bypass this problem 
has been proposed [9| and implemented experimentally ifTUl . 
The idea is to change the trapping potential in a way that 
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prepares the same final state as the adiabatic process at a given 
final time. This method, neatly characterized as "shortcut to 
adiabaticity" ifTTI . provides a family of paths which achieve the 
desired frictionless expansion, and in theory the necessary time 
can be made arbitrarily small. In practice, there are always 
experimental constraints which limit this time to a finite value. 
In this article, we impose restrictions motivated by exper- 
iments and formulate the desired transfer as a time-optimal 
problem. We then use optimal control theory to find the short- 
est adiabatic-like path for a one-dimensional BEC, trapped 
in a time-dependent harmonic potential. We show that the 
minimum time frictionless expansion takes place when the 
frequency of the potential changes in a bang-bang manner 
between the allowed values and calculate the necessary ex- 
pansion time. Note that numerical optimization methods have 
been used to control a BEC in an optical lattice [12] or in a 
magnetic microtrap |fi3l . Our approach here is different, since 
we use the time-optimal theory of single-input control systems 
in the plane [14), lfl5l . The analysis continuous our previous 
work |[T6l . IfTTI . where we have considered the minimum time 
frictionless cooling of a harmonically trapped atom. 

II. FORMULATION OF THE PROBLEM IN TERMS 
OF OPTIMAL CONTROL 

The evolution of the wavefunction tp(t,x) of a condensate 
trapped in a one-dimensional (elongated cigar trap) parabolic 
potential with time-varying frequency ui(t) is given by the 
following Gross-Pitaevskii equation [9| 
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where m is the particle mass, g is the coupling constant 
and h is Planck's constant; i e M and ip is a square- 
integrable function on the real line. When ui(i) is constant 
and g\ip\ 2 / (Hui) ^S> 1 the kinetic energy term may be neglected 
ifTHll . the so-called Thomas-Fermi approximation, 
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The above equation can be solved by separation of variables 
and the solution is 
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Fig. 1 . Time evolution of the harmonic trap frequency. 




and the chemical potential fi u is determined from the number 
of bosons N through the normalization condition 



\?p(t,x)\ 2 dx = / \y uJ (x)\ 2 dx = N. 



We find easily that 
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Consider now the case shown in Fig. [TJ where u)(t) = ljq 
for t < and ui(t) — lot < uq for t > T, For frictionless 
evolution, the path oj(t) between these two values should be 
chosen such that if 



then 



ip(fl,x) = * Uo (x), 



1>(t,x) = e ia ®* UT (x), t>T, 
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where a(i) is a global (independent of the spatial coordinate 
x) phase factor. According to © and ©, this evolution 
corresponds to the expansion of the condensate by the factor 
(ujo/ijJt) 2 ^ 3 , see Fig. 12 Among all the paths u(t) that result 
in (O, we would like to find one that achieves frictionless 
expansion in minimum time T. In the following we provide 
a sufficient condition on ui(t) for frictionless expansion and 
we use it to formulate the corresponding time-optimal control 
problem. 

Proposition 1: If ui(t), with w(0) = loq and ui(t) = lj(T) = 
ujt for t > T, is such that the equation 



b(t)+oj 2 (t)b(t) = 



b 2 (t) 
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has a solution b(t), with 6(0) = 1,6(0) = and b(t) = 
b(T) = (lu /uj t ) 2/3 , t>T, then condition ((7} for frictionless 
expansion is satisfied. 

Proof: The frequency variations in the trapping potential 
change the time and distance scales and motivate the use of 
the following "ansatz" fl9l, Ifl9ll, in fl) 
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where \ = x/b(t), t = r(t) is a time rescaling, and 
the distance scale bit) satisfies <[8J an d me accompanying 



Fig. 2. Schematic representation of the condensate expansion. With solid 
line we represent the initial harmonic potential (convex) and the corresponding 
wavefunction (concave), while with dashed line we depict the final potential 
(convex) and the expanded wavefunction (concave). 



boundary conditions. We obtain 
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with the initial condition <^>(0, %) — ^u ix)- l n tn i s equation 
the frequency is constant and we can apply the Thomas-Fermi 
approximation neglecting the kinetic energy term, provided 
that the scaling factor b does not become very small (we will 
secure this by imposing appropriate conditions later), 
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We will show that for t > T, where 6(t) = (w /wt) 2/3 , 
ipit,x) has the desired form. We examine separately each 
of the three terms in (fl2t . Since 6(i) = in this time 
interval, the first exponential is equal to unity. About the 
second exponential, observe from (TlOb that 

2/3 
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since 6(t) = (w /<^t) 2 ^ 3 , t > T. Also, from (O we have 
Mojo/Pwt = ( w oK) 2/3 - Thus 
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The last term in <TT~2T > satisfies 

1/3 



LOT 
UJq 



#., 



U>T 

LJQ 



2/3 



*wr(x), 



as it can be verified using (0|l and ©. Putting all these together 
we see that ip(t, x) has the desired form for t >T. ■ 

In order to find the path w(t), < t < T, that accomplishes 
frictionless expansion in minimum time T, we express the 
problem using the language of optimal control, incorporating 
possible restrictions on u(t) and b(t) due, for example, to 
experimental limitations. If we set 



xi = b, x 2 
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and rescale time according to tnew = ^o^old' we °b ta i n 
the following system of first order differential equations, 
equivalent to equation ([H} 
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If we set 7 = (uo/wr) 2 ^ 3 > 1, the time optimal problem 
takes the following form 

problem 1: Find -1 < u(t) < 1 with u(0) = l,u(T) = 
I/7 3 such that starting from (cci(O), #2(0)) = (1, 0), the above 
system reaches the final point (xi(T), x 2 (T)) — (7, 0), 7 > 1 
in minimum time T, under the constraint x 2 > 0. 

The boundary conditions on the state variables (xi, x 2 ) are 
equivalent to those for b, b, while the boundary conditions 
on the control variable u are equivalent to those for ui, so 
the requirements of Proposition Q] are satisfied. Note that the 
possibility ui 2 (t) < (expulsive parabolic potential) for some 
time intervals is permitted J5), iTTTI . The maximum allowed 
frequency is taken to be equal with the initial frequency 
uJq ||20l . The lower bound is taken as — lu 2 since the same 
harmonic potential can be made repelling or attractive by 
a phase shift when it is created by an optical lattice |20|. 
The assumption x 2 > implies that b > 0, so the scale 
factor is an increasing function of time satisfying b(t) > 1. 
This assumption guarantees that the condensate is continu- 
ously expanding and excludes breather-like solutions which 
might be difficult to implement experimentally. Note that 
the aforementioned condition holds for both time profiles 
of bit) used in |[9). Specifically, for the polynomial form 

b p( s ) = (7 - iX 6 ^ ~ 15s4 + 10s3 ) + 1 > where s = l l T - 

it is b p (s) = 30(7 - l)s 2 (s - l) 2 > 0, while for the 
exponential of polynomial form b e (s) = j 6s ' J - 15s + l0s it 
is b e (s) = 301n7s 2 (s - l) 2 b e (s) > 0. Now observe that the 
system (fT4l and (fl31 l describes the one-dimensional Newtonian 
motion of a unit-mass particle, with position coordinate xi and 
velocity x 2 . The acceleration (force) acting on the particle 
is — UX\ + l/x 2 . From this point of view, the breather- 
like solutions correspond to retrograde motion [21 1. The full 
benefits from this kind of motion arise when the particle 
moves close to the strong repulsive potential at X\ = 
|[T6l . But when b —> the Thomas-Fermi approximation in 
(fTTT > becomes questionable. For all these reasons we will not 



consider breather-like solutions here and we defer the study 
of retrograde motion to a future publication. 

In the next section we solve the following optimal control 
problem 

problem 2: Find —1 < u(t) < 1 such that starting from 
(xi(0),x 2 (0)) — (1,0), the system above reaches the final 
point (xi(T),x 2 (T)) = (7,0), 7 > 1 in minimum time T, 
under the constraint x 2 > 0. 

In both problems the class of admissible controls formally 
are Lebesgue measurable functions that take values in the 
control set [—1,1] almost everywhere. However, as we shall 
see, optimal controls are piecewise continuous, in fact bang- 
bang. The optimal control found for problem|2]is also optimal 
for problem [T] with the addition of instantaneous jumps at 
the initial and final points, so that the boundary conditions 
u(0) = 1 and u(T) = I/7 3 are satisfied. Note that in 
connection with Fig. Q] a natural way to think about these 
conditions is that u(t) = 1 for t < and u(t) — I/7 3 for 
t > T; in the interval (0, T) we pick the control that achieves 
the desired transfer in minimum time. This approach is similar 
to that used in our recently published work ifPTl . Il22l . as well 



III. OPTIMAL SOLUTION 
The system described by ( [Pil l. (fl3T l can be expressed in 



compact form as 

x = f(x) +ug(x), 
where the vector fields are given by 
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U = [—1,1]. Admissible controls are Lebesgue measurable 
functions that take values in the control set U. While the 
bound in x 2 is imposed to exclude retrograde motion, the 
bound in X\ results form the initial condition and the system 
dynamics. Specifically, observe that starting with any positive 
initial condition rri(O) > 0, and using any admissible control 
u, as x\ — > + , the "repulsive force" 1/x 2 leads to an increase 
in x\ that will keep X\ positive (as long as the solutions exist). 
Given an admissible control u defined over an interval [0, T], 
the solution x of the system ( [ToT l corresponding to the control 
u is called the corresponding trajectory and we call the pair 
(x, u) a controlled trajectory. 

For a constant Ao and a row vector A = (Ai, A2) € (M 2 ) 
define the control Hamiltonian as 

H = if (A , A, x, u) = A + (A, f(x) + ug(x)). 



Then, Pontryagin's Maximum Principle [24| provides the 
following necessary conditions for optimality: 

Theorem 1: Let (x*(i),Tt*(t)) be a time-optimal controlled 
trajectory that transfers the initial condition x(0) = xo into the 
terminal state x(T) = Xx. Then it is a necessary condition for 
optimality that there exists a constant Ao < and nonzero, 
absolutely continuous row vector function X(t) such that: 



1) A satisfies the so-called adjoint equation 

\(t) = - — (\ ,\(t),x*(t),u*(t)) 

= - (\(t),Df{x*(t)) + u,{t)Dg{x,{t))) (18) 

2) For < t < T the function u >->• H(\ ,\(t),x*(t),u) 
attains its maximum over the control set U at u = u^(t). 

3) jy(Ao,A(t),a?*(t),tt.(*))sO. 

We call a controlled trajectory (x, u) for which there exist 
multipliers Ao and X(t) such that these conditions are satisfied 
an extremal. Extremals for which Ao = are called abnormal. 
If Ao < 0, then without loss of generality we may rescale the 
A's and set Ao = —1. Such an extremal is called normal. 

Definition 1: We denote the vector fields corresponding to 
the constant controls u = — 1 and u = +1 by X = f — g and 
Y — f+g, respectively, and call the corresponding trajectories 
the bang X- and y-trajectories. A bang-bang trajectory is a 
finite concatenation of X- and F-trajectories. A concatenation 
of an X-trajectory followed by a F-trajectory is denoted by 
XY while the concatenation in the inverse order is denoted 
by YX. 

Proposition 2: For Problem[2]the optimal trajectory has the 
bang -bang form XY. 

Proof: For the system (TT4l i. (TT~5T > we have 



H(X , A, x, u) — A + \ix 2 + A 2 
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and thus 
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Observe that H is a linear function of the bounded control 
variable u. The coefficient at u in H is — A2X1 and, since 
x\ > 0, its sign is determined by $ = — A2, the so-called 
switching function. According to the maximum principle, point 
2 above, the optimal control is given by u = — 1 if $ < and 
by u = 1 if $ > 0. The maximum principle provides a priori 
no information about the control at times t when the switching 
function $ vanishes. However, if $(i) =0 and $(£) 7^ 0, then 
at time t the control switches between its boundary values 
and we call this a bang-bang switch. If $ were to vanish 
identically over some open time interval / the corresponding 
control is called singular. Now observe that, whenever the 
switching function $(£) = — A2(t) vanishes at some time t, 
then it follows from the non-triviality of the multiplier A(i) 
that its derivative <£>(£) = — Mit) — Xi(t) is non-zero. Hence 
the switching function changes sign and there is a bang-bang 
switch at time t. Thus optimal controls alternate between the 
boundary values u = — 1 and u = 1 of the control set. 

We next show how the restriction x 2 > limits the number 
and the type of switchings. For u = 1 the initial point (1, 0) 
is an equilibrium point for system (fl4] >. (fT5t so we should 
start with u = — 1, moving the system along an X-trajectory. 
Observe from (fT~5T > that for u = — 1 it is ±2 > so X2 > 0, 
and a switching to a y-trajectory at a point with x-i (t) > is 
necessary in order to reach the target point (7, 0). Then, since 
M(i) — at the switching point, it follows from H = that 
\i\t)x 2 {t) = -Aq. Since Ai(t) = $(<) > 0, it is A < 0. We 
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C(y,0) 



Fig. 3. One- and multi-switchings trajectories. Solid line corresponds to 
u = 1, dashed line to u = — 1. 



henceforth only consider normal trajectories and set Ao = — 1. 
Then, H = implies that for any switching time t' we must 
have X 1 (t')x 2 {t') = 1. Since x 2 (t') > it is $(£') = Ai(t') > 
and obviously the initial switching from X to Y is unique. 

■ 

An intuitive way to understand why the one-switching 
trajectory is the fastest is to compare it with a trajectory 
with more switchings, like that shown in Fig. [3] and use 
the particle model mentioned above. In both trajectories the 
particle traverses the same distance x\ but along the trajectory 
with one switching its speed x\ = x 2 is always higher. 

There is an elegant way to show the time-optimality of the 
one-switching trajectory, at least when compared to trajectories 
with multiple switchings at points with x 2 > 0, using the 
one-form tt introduced in [14| and defined by $!(/) = 1 
and £l(g) — 0. Note that for x 2 7^ the vectors f,g are 
linearly independent in D so the one-form can be defined in 
this two-dimensional manifold by its action on these fields. 
Now consider the parts where the two trajectories shown in 
Fig. [3] differ and denote them by c (part B\BB^) and c' 
(part BiB 2 B^BiB^). The infinitesimal displacement along 
each trajectory is dx = (f + ug)dt, so dt — il(dx). 
The necessary time to travel cist = J dt = J il and 
the corresponding time for c' is t' = J dt = J , fL So 
t' — t = r,0 — f fl = <f _, , fl, where cr 1 is c run 
backwards. Since c^ 1 * c' is oriented counterclockwise, by 
Stokes theorem we have t' — t = J R dCl, where R is the region 
enclosed by c" 1 * c'. From the definition of 51 we can find its 
expression in cartesian coordinates Q, = (l/x 2 )dxi and from 
this the exterior derivative dil — (l/x^dxi A dx 2 E5\ . Thus 

t' — t= / dtt = / —dxidx 2 > 
Jr J r x 2 

and obviously the trajectory with only one switching is faster. 
We now move to calculate the minimum time T necessary 
to reach the final point (7,0). A first integral of the motion 
along the X-trajectory starting from (1,0) is 

, 2 



x 2 
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(21) 



while a first integral of the motion along the F-trajectory 
ending in (7, 0) is 



The substitution 



x 2 +x 2 + — =7 2 + -, 
xi 7 



(22) 



The total time T is given by the sum of the times spent on 
each trajectory segment 



T = T 1 +T 2 . 



(23) 



Since x\ — X2 > for both cases, we can easily find from 



T-i = 
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(xi - l)(jcf + a;i +2) 



dx x (24) 



and from d22l 



T 2 = 



Xx 

p V (7- x-l){x\ +7*1-2/7) 



dxi, (25) 



where 5(/3, (5) is the common point of the X and F segments 
shown in Fig. [3] satisfying d2"TT i and (l22l . It is 
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We express the times Tx,T 2 in terms of elliptic integrals and 
elementary functions. We start from (l24l by rewriting the 
integrand in the convenient form (note that x\ 7^ 0) 



Ti = 
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dxx. (27) 



Following the procedure described in 11261 . we can bring the 
above integral in the form 
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allows us to express Ix in terms of elliptic integrals 
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are the incomplete elliptic integrals of the first and third kind, 
respectively, in Jacobi's form, while 
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On the other hand, the substitution 
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gives 
h = 



c=(l+p 2 )(l-q 2 ), 



dw 
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2^7i V^^T 2 \/c 
Using all the above, the time 7\ spent on the X segment of 
the trajectory can be calculated as a function of 7. We move 
on to calculate T%. Integral d25l l can be written as 
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where 



7 + V7 2 + 8/7 



c = z2±VjHNl. (45) 



Since — e < < £ < 7 and £ < /3 < 7 (recall that 7 > 1), we 
find from E3 that 
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The integral in d28l i can be decomposed as 
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V(7-O(0 + e) 
Note that /3 is given in d26i > as a function of 7. In Fig. [4] we 
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plot the total expansion time T = Tx + T 2 with respect to 7. 
We finally show that for large 7 time T grows logarith- 
mically, which seems to be a universal characteristic for 
processes involving shortcuts to adiabaticity ifTBI , ||20) , ||28l . 
Instead of taking the limit 7 — >• 00 in the above expressions, 




r 



Fig. 4. Total time for the adiabatic-like expansion as a function of the 
expansion factor. 



we follow an easier path. Note that for large x\ the system 
equations ( [T4l and (TT3T > take the simple form 



Xl = X2, 
±2 = —UX\. 



(49) 
(50) 



For u = — 1 this system corresponds to a harmonic oscillator 
with angular frequency w = 1. The switching point i? in Fig. 
[3] tends to B -> (7/V2, 7/V2), so it is close to the line 
through the origin which makes a 7r/4 angle with the xi- 
axis. The necessary time to reach the final point (7, 0) from 
B approaches T2 — > 7r/(4w) = 7r/4. For u = 1 we find from 
(|49l and < T50b that for large xi and £ it is x\ ~ e*, so Ti ~ In 7. 
Thus, the dependence of total time for large 7 is T ~ In 7. 



IV. CONCLUSION AND FUTURE WORK 

In this paper, we combined optimal control theory with 
shortcut to adiabaticity to derive the optimal time variation of a 
confining harmonic potential which results in a fast adiabatic- 
like expansion of a trapped one-dimensional BEC. This fast 
frictionless expansion is an important task for applications 
where controlling BEC in an efficient way is crucial, for 
example atom interferometry. In our future work we would 
like to relax some of the assumptions used in this article, for 
example to permit breather-like solutions for the condensate 
evolution. 
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